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We show that chiral symmetry can be broken spontaneously in one-component systems with 
isotropic interactions, i.e. many-particle systems having maximal a priori symmetry. This is achieved 
by designing isotropic potentials that lead to self-assembly of chiral surfaces. We demonstrate the 
principle on a simple chiral lattice and on a more complex lattice with chiral super-cells. In addition 
we show that the complex lattice has interesting melting behavior with multiple morphologically 
distinct phases that we argue can be qualitatively predicted from the design of the interaction. 
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Breaking of chiral symmetry plays a central role both 
in fundamental physics, e.g., parity violation in the 
weak interaction, and in biology where many types of 
biomolecules exist only as enantiomers. The commonly 
accepted explanation for homochirality in chemical and 
biological systems is spontaneous symmetry breaking. 
This phenomenon involves two steps: first the formation 
of chiral molecules, and then a chiral specific catalysis 
that amplifies a stochastic imbalance to a macroscopic 
scale. In heterogeneous systems it is not hard to imagine 
that both steps can be achieved, and there are indeed 
many chemical systems that spontaneously deviate from 
racemic mixture []], Q- There are even several known 
molecules that show chiral auto-catalytic activity 0, Q . 
It is an interesting question whether spontaneous break- 
ing of chiral symmetry can happen also in simpler sys- 
tems, by which we here mean systems with higher degree 
of symmetry. The perhaps simplest such class of systems 
would consist of a single particle type with an isotropic 
pairwise interaction. For several reasons it appears to 
be much harder to form homochiral states in such sys- 
tems. First, cluster formation, which mimics the creation 
of molecules, typically results in achiral structures with 
dihedral group symmetries |5j. Second, (auto-) catalysis 
is not easily achieved in homogeneous isotropic systems. 
However, chiral symmetry can also be broken in another 
context, namely during crystallization. Again there are 
many examples of heterogeneous systems where this is 
known to happen @, Q, but no chiral crystals arising 
from isotropic interactions have been reported, neither 
in explorative [8| nor design [9| studies. In this Letter 
we show that one-component systems with carefully de- 
signed isotropic interactions can self-assemble into chi- 
ral lattices. Our results demonstrate that spontaneous 
breaking of chiral symmetry can occur in many-particle 
systems with maximal degree of symmetry. 

Chiral surfaces are of practical importance because of 
their potential for use in chiral catalysis. The most im- 
portant applications are in the pharmaceutical industry 
where we learned the hard way that a therapeutically 



well behaved enantiomer may be toxic in the other enan- 
tiomeric form. The classic example is the tragedy with 
Thalidomide, whose racemic form turned out to cause 
birth defects after being extensively used as a sedative for 
pregnant women [10]. Today the dominating method for 
synthesizing large quantities of chiral products involves 
various types of chiral catalysts in solution (homogeneous 
catalysis) [jj], fl2| - However, the development of chiral 
surface catalysts (heterogeneous catalysis) receive much 
interest due to their separability and reusability. Sev- 
eral ways of producing chiral surfaces exist, for example 
through surface reconstruction induced by adsorbed chi- 
ral molecules [13], cleaving achiral crystals along planes 
of lo w sy mmetry such as the (643) surface of an fee struc- 
ture [14], or self-assembl y o f molecules into chiral struc- 
tures on achiral surfaces One of the outstanding 
challenges for all heterogeneous chiral catalysis is how to 
make chiral surfaces with sufficiently large (macroscopic) 
active areas [17j . Simple models that exhibit formation of 
chiral domains, such as the models presented here, could 
prove useful in supporting progress in this area. 

In this work we focus on self-assembly at low tempera- 
ture where structure formation is driven by minimization 
of the potential energy. For a given configuration in a sys- 
tem with isotropic interactions, the energy is defined by 
the density distribution p(r) and is naturally described as 
a quadratic form that sums all the contributions from the 
pairwise potential. For our purposes it is suitable to ex- 
press the energy in reciprocal space where the quadratic 
form is diagonalized due to the translational invariance 
of the isotropic interactions [18| : 

E= y*dridr 2 p(ri)y(|n-r 2 |)p(r 2 ) 
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where \p(k)\ 2 dk = Jj sk | =fc elk |p(k)| 2 , p(k) is the stan- 
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FIG. 1: (a) The simplest form of chiral lattice, composed of 
scalene (acute) triangles with well separated side lengths U. 
The emerging chirality can be either left or right oriented (b) . 
(c) By selecting for the reciprocal lattice (the red peaks) of the 
target structure, a potential (d) which causes self-assembly 
into one of the two possible chiralities is obtained. 



Jrdr V(r) Jo(rfc). The small \k\ region describes the 
defining features of the crystal structure and by designing 
the energy spectrum of this region particles can be made 
to self-assemble into crystalline structures at low tem- 
peratures. We used this observation in a recent study 
to show how to design the interactions of a system so 
that it self-assembles into target lattices [20j |. Briefly, the 
method works as follows. The symmetry of the lattice 
manifests in reciprocal space as a restriction of the sup- 
port of p(k) to a set of finite points {G^}. By choosing 
V(k) smooth and positive with zeros coinciding with the 
reciprocal lattice of the target structure V(|Gi|) = 
[see Figs.QJc) and 12(c)], we can guarantee the target con- 
figuration to be a ground state. This construction utilizes 
the fact that there are a finite number of Bravais lattices, 
all with different structure factors. The structure factors, 
p(G^), will not affect the energy and hence all crystals 
with the same periodicity will be ground states. 

Chirality on a crystal surface is by definition equiva- 
lent to the absence of axes of symmetry and it can man- 
ifest itself in fundamentally different ways. Every ideal 
crystal structure consists of a mathematical lattice and 
a basis. [24] If either the basis or the lattice is chiral, 
the crystal will be as well. But even if both lattice and 
basis are achiral, the crystal can be chiral if the axes 
of symmetry of the lattice does not coincide with those 
of the basis. In this case the reciprocal lattice repre- 
sentation is also achiral and instead the chirality, if any, 
is determined by the distribution of the structure factors 




FIG. 2: Snub hexagonal tiling, a chiral Archimedean tiling. 
Note that the chirality does not originate from the positions 
of the hexagons (they are organized in rombs and not paral- 
lelograms as they must when the tiles are regular polygons) 
but from the orientations of the triangles in between, (a-d) 
as in Fig. [1] In (a) two regions of opposite chirality are high- 
lighted, (c) illustrates the perturbation of magnitude e at the 
maximum of p. 



(weights) over the reciprocal lattice points. We exemplify 
these fundamentally different chiralities with two lattices 
(crystal structures), a lattice of scalene triangles and the 
snub hexagonal tiling. For each case we discuss how the 
chirality relates to the reciprocal dual of the crystal and 
how this reflects on the energy spectrum of the poten- 
tial designed to self-assemble into the target structure. 
We also demonstrate the crystallization with snapshots 
of low energy states self-assembled when simulated in the 
canonical ensemble. 

The simplest chiral geometric form in two dimensions 
is the scalene triangle, where the chirality depends on the 
clockwise order of the side lengths. Hence the simplest 
chiral crystal is one forming scalene (acute) triangles, an 
oblique Bravais lattice. The oblique lattice is unique in 
the sense that it is the only chiral Bravais lattice. In this 
case the chirality of the lattice is directly manifest in the 
reciprocal lattice and will form scalene triangles with the 
same chirality as the reciprocal. Since the basis is trivial 
there is no freedom in the arrangement of its constituent 
particles and the only requirement for self-assembly of 
such a chiral structure is that the positive spectrum V(k) 
selects for the reciprocal lattice of the target and that the 
system has approximately the correct particle density. 
In Fig. [1(c) a spectrum fulfilling these criteria is shown, 
together with the corresponding potential and the result 
of a Monte Carlo simulation with the potential. We see 
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that the mirror symmetry of the system is indeed broken 
with a homochiral lattice as the result. 

As an example of a crystal where the chirality instead 
emerges from the interplay between the basis and the 
lattice, we use the snub hexagonal tiling and the crystal 
formed by its vertices, shown in Fig. [2j The tiling to- 
gether with its dual are the only uniform chiral tilings 
of the Euclidean plane. This pattern closely resembles 
geometries observed in experimental systems involving 
(anisotropic) molecules with a triangular geometry [2lj . 
An achiral superstructure of vacancies determines the lat- 
tice periodicity while the orientation of the basis of six 
particles in relation to adjacent vacancies determines the 
chirality. Unlike the previous example, here perturba- 
tions of the energy spectrum at |G;| are necessary to 
distinguish between crystals with different bases, since 
they consists of more than one particle; see [20] for more 
details on use of perturbations to break degeneracies of 
the ground state. There are six pairs of reciprocal lattice 
points representing wave vectors of length |G;| = \/7|ko|, 
ko being a primitive reciprocal vector. The Fourier trans- 
form of the crystal p(k) at these sites is equal to —1 or 
pmax = 6. The chirality is determined by the order of 
the structure factors in the pairwise sites. A negative 
perturbation — e of the energy spectrum at this k will 
ensure that the ground state will have the correct basis, 
as shown in Fig. Efc). In the zero temperature limit, a 
small deformation in the form of a small rotation of the 
basis by 0.2 rad towards its chiral counterpart will be 
visible: the spectrum will not be affected at its maxima 
by small perturbations of the basis, while a rotation in- 
creases the negative energy contribution from the other 
peaks at that wavelength. The maxima correspond to 
the hexagonal close-packed lattice which, together with 
the vacancies introduced by the limited density, forms 
the snub hexagonal tiling. 

To test whether the target lattice can be assembled 
from random initial configurations we perform Monte 
Carlo simulations of particles interacting with the ob- 
tained potentials, allowing only local moves and with 
simulated annealing at fixed density. During the tem- 
perature annealing the particles organize into homochi- 
ral regions. The boundaries where the chirality changes 
are energetically disfavored and the homochiral grains 
grow as the system is annealed towards one of its ground 
states. The annealing process is however very slow, as 
illustrated in Fig. [3] where the time evolution is shown 
on an (approximately) logarithmic time scale. Slow do- 
main growth is typically observed also in experimental 
systems, which is causing practical difficulties in the syn- 



thesis of chiral surfaces [17 ]. 



Further simulations reveal that the snub hexagonal 
model has a rich phase diagram, as demonstrated in 
Fig. |U It shows a large region of stability for the target 
chiral configuration (c) as temperature and magnitude of 
the spectrum perturbation e (see Fig. [2]) is varied. We 



can understand the various neighboring phases in terms 
of the target structure losing some of its properties while 
retaining others. If the negative perturbation stabilizing 
the basic hexagonal lattice is large, the preference for par- 
ticles to occupy sites on this lattice will remain even at 
temperatures where the periodicity of the snub hexago- 
nal lattice, induced by the locations of zeros in the energy 
spectrum, disappears. The crystal melts through a va- 
cancy unbinding resulting in a hexagonal lattice with ran- 
domly placed vacancies, Fig. HJc). For even larger e the 
vacancies aggregate into clusters, Fig. ffljd), not forming 
the snub hexagonal lattice for any /}. If the negative per- 
turbation is small, the basis will dissolve at temperatures 
below those where the periodicity of the lattice is broken. 
A large variety of basis configurations are observed in this 
region although the dominant and most stable morphol- 
ogy is striped, Fig. Sfb). We suspect that the rich phase 
behavior observed here is typical for models constructed 
through our design method as the spectrum fixes a set of 
features of the ground state with varying strengths, each 
breaking down under different conditions. This demon- 
strates that the energy spectrum approach to designing 
potentials for targeted self-assembly can also be useful 
when trying to understand and control complex finite 
temperature phase behavior. 

In summary, we demonstrate that isotropic pairwise 
potentials can cause particle systems to self-assemble into 
chiral surfaces with varying degree of complexity. In gen- 
eral this work hints at some of the future possibilities of 
nano scale self-assembly, expanding on what we know to 
be possible, and could be an inspiration for the material 
science of tomorrow. Lately, there has been remarkable 
progress on designing nano and colloidal particles with 



exotic interaction potentials 22, 23]. Still, experimental 



realizations of systems with as complicated interactions 
as we use here are not likely to appear in the near future. 
For theoretical methods and experimental techniques to 
converge and the field of self-assembly moving from sci- 
ence to industry on a larger scale, further progress from 
both directions is necessary. 

OL and MNJ acknowledge support from the SuMo Bio- 
materials center of excellence. 



* Electronic address: mjacobi@chalmers.se 
[1] Y. Okamoto and T. Nakano, Chem. Rev. 94, 349372 
(1994). 

[2] E. Eliel and S. Wilen, Stereochemistry of Organic Com- 
pounds (Wileyi-Interscience, 1994). 

[3] K. Soai, T. Shibata, H. Morioka, and K. Choji, Nature 
378, 767 (1995). 

[4] K. Soai and I. Sato, Chirality 14, 548 (2002). 

[5] E. Edlund, O. Lindgren, and M. N. Jacobi, Phys. Rev. 
Lett. 107, 085501 (2011). 

[6] D. K. Kondepudi, R. Kaufman, and N. Singh, Science 





e 




FIG. 4: A phase diagram for the snub hexagonal model shows 
the equilibrium configuration as a function of inverse temper- 
ature /3 and the magnitude e of the perturbation. Without 
the perturbation, striped and stripe-like morphologies are ob- 
served (a), a strong perturbation leads to vacancy unbinding 
(b) before complete disorder (c) as temperature increases. A 
too strong perturbation results in aggregation of the vacan- 
cies (d). In the central region (e) the chiral snub hexagonal 
lattice is formed. 



250, 975 (1990). 
[7] H. Flack, Helv. Chim. Acta 86, 905 (2003). 
[8] M. Engel and H.-R. Trebin, Phys. Rev. Lett. 98, 225505 

(2007). 

[9] S. Torquato, Soft Matter 5, 1157 (2009). 



[10] G. Blaschke, H. Kraft, K. Fickentscher, and F. Kohler, 
Arznei.-Forschung 29, 1640 (1979). 

[11] T. P. Yoon and E. N. Jacobsen, Science 299, 1691 (2003). 

[12] E. N. Jacobsen, A. Pfaltz, and H. Yamamoto, Compre- 
hensive Asymmetric Catalysis, Vol. 1-3 (Springer, Berlin, 
1999). 

[13] M. Schunack, L. Petersen, A. Kiihnle, E. Laegsgaard, 

I. Stensgaard, I. Johannsen, and F. Besenbacher, Phys. 

Rev. Lett. 86, 456 (2001). 
[14] D. S. Sholl, A. Asthagiri, and T. D. Power, J. Phys. 

Chem. B 105, 4771 (2001). 
[15] M. O. Lorenzo, C. J. Baddeley, C. Muryn, and R. Raval, 

Nature 404, 376 (2000). 
[16] H. Spillmann, A. Dmitriev, N. Lin, P. Messina, J. V. 

Barth, and K. Kern, J. Am. Chem. Soc. 125, 10725 

(2003). 

[17] A. J. Gellman, ACS Nano 4, 5 (2010). 
[18] E. Edlund and M. N. Jacobi, Phys. Rev. Lett. 105, 
137203 (2010). 

[19] G. B. Folland, "Fourier Analysis and Its Applications," 
(Brooks/Cole Pub Co, Pacific Grove, 1992) Chap. 7.5, p. 
246. 

[20] E. Edlund, O. Lindgren, and M. N. Jacobi, Phys. Rev. 
Lett. 107, 085503 (2011). 

[21] S. Katano, Y. Kim, H. Matsubara, T. Kitagawa, and 
M. Kawai, J. Am. Chem. Soc. 129, 2511 (2007). 

[22] C. N. Likos, Phys. Rep. 348, 267 (2001). 

[23] Y. Min, M. Akbulut, K. Kristiansen, Y. Golan, and 
J. Israelachvili, Nat. Mater. 7, 527 (2008). 

[24] We have up to now followed the common but inexact 
usage of the word lattice to mean any ideal crystal struc- 
ture, but in the following three paragraphs we need to be 
more exact. 



